Thermodynamic equivalence of certain ideal Bose and Fermi gases 
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We show that the recently discovered thermodynamic "equivalence" between noninteracting Bose 
and Fermi gases in two dimensions, and between one-dimensional Bose and Fermi systems with 
linear dispersion, both in the grand-canonical ensemble, are special cases of a larger class of equiv- 
alences of noninteracting systems having an energy-independent single-particle density of states. 
We also conjecture that the same equivalence will hold in the grand-canonical ensemble for any 
noninteracting quantum gas with a discrete ladder-type spectrum whenever aA/NksT is small, 
where N is the average particle number and a its standard deviation, A is the level spacing, fce is 
Boltzmann's constant, and T is the temperature. 

PACS numbers: 05.30.-d, 05.70.Ce 
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I. INTRODUCTION 

There has been considerable recent interest in surpris- 
ing thermodynamic "equivalences" between certain ideal 
Bose and spinless Fermi gas systems, including nonrcl- 
ativistic free particles in two dimensions^ and between 
one-dimensional particles with a linear, sound-like dis- 
persion relation^ Both results are valid in the grand- 
canonical ensemble, and assert that the Hclmholtz free 
energies F F and F B of the Fermi and Bose systems, re- 
spectively, are simply related by 

N 2 
F F (T,V,N)-F B (T,V,N) = —, (1) 

where C is a constant that may depend on the system 
volume V, but is independent of temperature T and the 
mean particle number N. The Bose and Fermi systems 
are assumed to have identical single-particle Hamiltoni- 
ans, and have the same T, V, and N. 

There are several immediate consequences of Eq. {Q, 
including: 

a. The entropies of the Fermi and Bose systems are 
identical, 



S F {T,V,N)=S B (T,V,N); 



(2) 



b. Their internal energies differ by a temperature- 
independent constant, namely 

N 2 
U F (T, V, N) - 17 B (T, V,N) = —; (3) 

c. The constant-volume heat capacities are identical; 

d. The chemical potentials are simply shifted by a 
temperature-independent constant, 



e. The thermodynamic potentials are connected by a 
relation opposite to JD, 



N 2 
n F (T,V,N)-Q B (T,V,N) = - — ; (5) 



f. The pressures of the Fermi and Bose gases satisfy 

N 2 



P F (T,V,N)-P B (T,V,N) 



2CV 1 



(6) 



and again differ only by a temperature-independent 
constant. 

These results explain and considerably extend isolated 
thermodynamic relations that were discovered some time 
ago by May£ 

In this paper we demonstrate that the equivalence de- 
fined in Eq. Q holds whenever the single-particle den- 
sity of states (DOS) is independ ent of energy, for which 
the systems considered in Refs. Ill2u .il are special cases. 
We then extend recent results of Schmidt and Schnack^ 
and of Crescimanno and Landsberg^ on harmonically 
confined Bose and Fermi gases in one dimension to fur- 
ther enlarge the class of equivalences to include any ideal 
quantum gas with a discrete ladder-type spectrum, in the 
limit 



gA 
Nk B T 



0. 



(7) 



TV 
fi F (T,V,N) - f, B (T,V,N) = -; 



(4) 



Here iV is the mean particle number, a is the standard 
deviation in particle number about N, A is the energy- 
level spacing, k B is Boltzmann's constant, and T is the 
temperature. The condition Q requires that either (i) 
the particle number is conserved, (ii) the number of par- 
ticles is very large, (iii) the temperature is very high, (iv) 
the spectrum is continuous, or any combination of these 
possibilities. 



II. THERMODYNAMIC EQUIVALENCE AND 
THE DENSITY OF STATES 



We begin by writing the grand-canonical partition 
function of an arbitrary nonintcracting Bosc or Fermi 
system as 



a N a 



(8) 



Here a labels the quantum states of a single Bosc or spin- 
less Fermi particle with spectrum e a , and j3 = l/fc B T. 
The occupation numbers N a take the values N a = 
0, 1, 2, • • • for bosons and N a = 0, 1 for fcrmions. The 
thermodynamic potential 51 = F — fj,N is given by 



fi= —=ha.Z. 



(9) 



Because the average number of particles is required to 
be the same for the Bose and Fermi cases, their chem- 
ical potentials \i in Eq. (JSJ arc different. The relations 
between // B , MF; and N, are determined by 

^n B (e Q -Mb) = ^n F (e Q - /j f ) = N, (10) 



where 

n B (x) 



and n F (x) 



; /3x _ i 



e^ + 1 



(11) 



are the Bose and Fermi distribution functions. 
Next we define a single-particle DOS according to 



above that minimum. Without loss of generality we can 
take the minimum to be at e = 0. Then 



ff(e)=Ce(c), 



(16) 



where 0(e) is the unit step function. The most common 
example of a DOS of the form (|16[) occurs for free nonrel- 
ativistic particles of mass m in two dimensions, in which 
case 

2771 
C =^2- 17 

However, there are other situations where (|16fl holds as 
well, including nonintcracting particles moving in one di- 
mension with a linear dispersion e(k) oc |fc|, and also 
for particles moving in three dimensions with qubic dis- 
persion e(fc) ex \k\ 3 . These cases were noted earlier by 
PathriaiSi Furthermore, we note that the equivalence (|TJ 
would not apply to two-dimensional systems moving in 
the potential of a corrugated surface or to ideal lattice 
gas models. 

Assuming (116(1 . we can immediately obtain 



C_ 

~ 2 



^B = -^Li 2 {zb) 



and 



C_ 

~ 2 



^f = : ^Li 2 (-z F ), 



(18) 



(19) 



where z B = 1 - e ~ pN / c and z F = e fm ' c - 1. Here 
\A<z(x) is Euler's dilogarithm function. Furthermore, 
from Eq. (|15|) . we have /i F — fJ-B = N/C and z-p = 
zb/(1 — zb)- These relations, along with the identity 



3(e) = X^( £ ~ ec ")' 



(12) 



which gives the number of energy levels per unit energy, 
as a function of e. In a translationally invariant system, 
g(e) scales linearly with system volume V. In terms of 
the DOS, wc have 

1 f°° 

^B = ~3 deg(e)ln(l-e-^ZB) (13) 

P J — oo 

and 

1 r°° 
Qf = —s <kg(e)\n.(l + e- f,e z F ), (14) 

P J —oo 

where zb = eP^ and z F = e^^ F are the Bose and Fermi 
fugacities. Furthermore, condition l|lUf) can be written as 



de 



?(e) 



e^z^-1 



de 



9(e) 



e^ zp l + 1 ' 



(15) 



The expressions <|13ll . I|14|) . and i|15|) . are valid for any 
nonintcracting quantum gas. 

We assume that spectrum is bounded from below, and 
that the DOS is a constant, C, independent of energy, 



Li 2 ( a; )+Li 2 (^ T ) = -i[ln(l-x)] i 



(20) 



directly lead to the equivalence stated in Q. The ther- 
modynamic equivalence evidently applies to any nonin- 
tcracting quantum gas with a constant DOS. 

It is also instructive to directly demonstrate the equiv- 
alence of the entropies: For a system with a constant 
DOS of the form (|16(1 . the Bose and Fermi entropies ara& 



Sb = —Ckv 



de 



n B (e-/X B ) ln[n B (e-/j B )] 



[l + n. B (e-/i B )] ln[l + n B (e-/i B ) 



(21) 



and 



5 1. 



-Cfc B / de 
Jo 



n F (e-/j F )ln[n F (e-/j F )] 



[l - n F (e-// F )] ln[l - n F (e-/j F ) 



(22) 



Changing the integration variable in the Bose case to 



w = e 13 ^ MB - > — 1, and in the Fermi case to w 
leads to 



Sb = CklT 



dw 



ln(l + w) 



In 



w 



o 0(z-Hf) 



(23) 
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FIG. 1: Helmholtz free energy per particle for a ID ideal Bose 
gas in a harmonic confining potential with level spacing A. 



and 



S F = Ck&T 



dw 



ln(l 



lnu> 
1 +w 



(24) 



Notice that the statistics dependence enters only in the 
lower integration limits. Because the average particle 
numbers are the same, these lower limits coincide and 
thus the Fermi and Bose entropies are identical. 



III. ID QUANTUM GASES IN HARMONIC 

POTENTIALS AND AN EXTENDED 

EQUIVALENCE 

It is interesting to consider whether a constant DOS is 
necessary for the equivalence defined in Eq. Q . In partic- 
ular, does it apply to a system with a discrete ladder-type 
spectrum of the form 



e n = nA, n = 0,l,2,3,. 



(25) 



where A is the level spacing, which reduces to the case 
considered in Sec. [n] in the limit A — > 0? In this section 
we demonstrate that the equivalence does still hold, in 
the grand-canonical ensemble, for noninteracting gases 
with the spectrum l|25|) . in the limit N — > oo, and also in 
the canonical ensemble for any N. 

In Figs, ^and |21 we show the Helmholtz free energy per 
particle, numerically calculated in the grand-canonical 
ensemble, for ID quantum gases in a harmonic potential 
with level spacing A. The free energies and temperatures 
are plotted in units of A. In these figures, the solid curves 
are for TV = 1000 and the dashed curves are for N = 
10. In Fig. |21 the difference between the Fermi and Bose 
free energies are given as a function of temperature, with 
Fp shifted by the Fermi ground-state energy E N [defined 
below in Eq. i|30|l ] for convenience. 

For small N (dashed curve in Fig. I3J , the free energies 
clearly do not differ by a temperature-independent con- 
stant. However, as N becomes larger (solid curve), the 
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FIG. 2: Helmholtz free energy per particle for a ID ideal 
Fermi gas in a harmonic confining potential with level spacing 

A. 



equivalence does apply. These numerical results suggest 
that when particle-number fluctuations become negligi- 
ble, the equivalence holds. To establish this result, we 
use the fact that in the large- N limit, the grand-canonical 
free energy approaches the canonical free energy (shown 
in Figs. H an dEl as dotted curves), and that the equiva- 
lence holds exactly in the canonical ensemble. 

Writing the grand-canonical partition function Z in 
terms of the canonical partition functions Z N according 
to 



Z = J2 Z ^, * 



= e^\ 



N=l 



where z is the fugacity, leads to 

1 fd N Z 



N\ \ dz 



(26) 



(27) 



For Bose (C = l) and Fermi (C = — 1) particles with spec- 
trum H25[l . the grand-canonical partition function is 



Z = cxp 



-cxxi-cy* 



with b = e~ l3A , (28) 



from which we obtain 



Z N = e~ 0EN x 



ife) 



(29) 



Here 



EL 





N(N-l) 



A 



for bosons 
for fermions 



(30) 



is the ground-state energy of N particles. The result 
in Eq. I|29|) was also obtained by Schmidt and Schnack 4 
using related methods. 
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FIG. 3: Free energy difference per particle for ideal quan- 
tum gases in a harmonic confining potential. Here E N is the 
ground-state energy of the Fermi system. Not shown is the 
difference in canonical ensemble case, which is exactly zero 
for any N. 



According to Eq. i|29|) , the canonical free energy for N 
particles is simply 



N 



F N = £ N + fc B T^ln(l 



-pAj-^ 



(31) 



scmblc is the exact mapping between the many-particle 
excitation spectra of both systems. 

Finally, we note that the thermodynamic equivalence 
with the spectrum (|25J) also trivially holds in the T — > oo 
limit, because in this limit the systems become classical. 
Therefore we conjecture that the equivalence defined in 
Eq. £Q| will hold in the grand-canonical ensemble for any 
nonintcracting quantum gas with a discrete ladder-type 
spectrum, whenever the quantity 



aA 

Nk B T 



(32) 



is small, where N is the average particle number and a is 
its standard deviation about N. This ratio roughly char- 
acterizes the magnitude of energy fluctuations caused 
by the exchange of particles with the environment — if 
allowed — relative to the thermal energy. We conclude 
that the (partial) thermodynamic equivalence discovered 
by Leei holds for ideal quantum gases with the spectrum 
(|25[) whenever (i) the number of particles is strictly con- 
served, (ii) the number of particles becomes very large so 
that a/N — > 0, (iii) T — > oo, (iv) A — > 0, or when any 
combination of these criteria are fulfilled. 
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